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ABSTRACT 


Photoelasticity is an optical method of experimental stress analysis which yields a whole- 
field representation of the principal-stress difference and its orientations. Though a large 
class of problems can be solved using the information of difference in principal stress alone, 
a certain class of problems need the value of individual stress components. The process of 
obtaining the individual values of the principal stresses separately is referred to as stress 
separation. Several separation methods have been developed to determine the individual 
values of the principal stresses. 

The present work discusses two methods of stress separation namely: shear 
difference method and Trebuna’s method for stress separation. Before the advent of 
widespread use of digital computers, shear difference technique was widely used for stress 
separation. The purpose of the present thesis is to have a relook on the applicability of this 
technique in digital photoelasticity. The various issues related to the digital implementation 
of the techniques are discussed in this thesis. To verify the implementation, initially fringe 
order and principal stress direction for the problem of a disc under diametral compression 
are used from theory of elasticity solution. The results compare very well with theory. 
Experimental data for the problem of disc under diametral compression is obtained using 
the phase shifting technique. As the phase shifting technique generally gives better values 
of isochromatic fringe order than the isoclinic value, initially stress components are obtained 
using experimental isochromatic fringe order and theoretical isoclinic value. Quantitative 
results for few arbitrary lines are compared with theoretical solutions. To verify the whole 
field implementation pseudo fringe patterns were plotted for individual stress components 
and compared with theory. 

In order to fully utilize the experimental data, the suitability of various methods for 
evaluating isoclinic value by phase shifting technique is discussed. The use of multiple load 
method for isoclinic determination, provided better values of isoclinic parameter for stress 
separation. Separation of stress components is then carried out for the problem of a finite 
plate with a hole in uniaxial tension. 
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Chapter 1 


INTRODUCTION 


1.1 Introduction 

Photoelasticity deals with determination of stress fields by the use of polarised light. 
This method is based on the temporary or artificial birefringence effect. Certain 
transparent, non crystalline materials like Plexiglas, celluloid, homolite, epoxy 1 etc., are 
optically isotropic when not loaded but become optically anisotropic and display 
characteristics similar to crystals when they are stressed 2 . The birefringence in the 
material is retained only during the application of the loads and disappears when they 
are removed. Measurements of this birefringence provides data adequate for quantitative 
determination of the corresponding state of stress. 

The measurements are made with the help of an optical instrument called the 
polariscope. A polariscope could be used to find out, namely, the maximum shear stress 
i.e. the difference in principal stresses and the orientation of principal stresses. The 
corresponding contours seen by the polariscope are called isochromatics and isoclinic 
fringes. The experimentation in photoelasticity involves acquisition of these images and 
their interpretation. To interpret the photoelastic images, data should be collected from 
the images and these data should be analyzed. The main merit of photoelasticity is its 
simplicity in experimentation as well as in the evaluation of the stress field. 
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1.2 Present work 

It is well known that a direct photoelastic approach cannot solve completely a general 
stress problem. The conventional photoelastic method provides two independent 
equations for a two-dimensional problem and five independent equations for a three 
dimensional problem 3 . Since there are three unknown stress components for a two 
dimensional problem and six unknown stress components for a three dimensional 
problem, the equations obtained by photoelasticity are not sufficient to determine the 
complete state of stress at any interior point. Hence, with a view to obtain additional 
relation, several auxiliary methods have been developed both for two and three 
dimensional problems, which in conjunction with photoelastic analysis provides the 
complete information for solving a general stress problem. Since all such methods 
provide additional information and enable the evaluation of stress components 
separately, these are referred to as stress separation methods. Each of these methods 
have their own advantage, limitation and application. 

The present work focuses on the shear difference technique. This method is 
relatively easier to employ for both two and three dimensional problems. This also 
happens to be one of the most widely used method, however it inherently bears the 
property of accumulation of errors, which further go on increasing, as one proceeds 
along the line of interest. Trebuna has modified the shear difference technique to 
minimize error accumulation. This method works on the principle of Block Integration. 
The isoclinic parameter which is the chief source of error is averaged along the 3 lines 
of the grid (the grid being employed as an integral feature of these methods) and 
assumed to represent the isoclinic along the centre line. With the noisy points of 
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isoclinic data being "smoothened out", this method is supposed to offer better results. 

A comparative analysis of these two methods is done for a disc under diametral 
compression. For this purpose first an error free set of data is theoretically generated 
through the standard solution of the stress field of circular disc under diametral 
compression. Trebuna’s method hardly shows any improvement over the conventional 
shear difference technique Error is intentionally introduced in the value of isoclinic 
parameter and it is observed that the averaging effect in Trebuna’s method does improve 
the accuracy of evaluation. 

Then use of experimentally obtained photoelastic data is studied. Use of phase- 
shifting technique, which uses the intensity data in spatial domain to extract the 
isochromatics over the field, is employed. In this technique, one has to record a few 
images corresponding to different optical arrangements for a given experimental set up. 
Phase-shifting technique, in general provides the estimation of fringe orders with 
reasonable accuracy. However the isoclinic parameter is not quite accurate. Hence, 
different optical arrangements to evaluate the isoclinic parameter such as a circular 
polariscope and a plane polariscope arrangement, are studied. The isoclinic parameter 
obtained by plane polariscope arrangement is found to be more accurate as against the 
circular polariscope arrangement and hence is suggested for analysis. The analysis is 
further improved by using multiple loads in plane polariscope arrangement for obtaining 
isoclinic data. 

With the two methods established for the disc under diametral compression, in 
the light of various combinations of photoelastic data, a similar analysis is performed 
for a plate with a hole, subjected to uniaxial tension, using completely experimental 


data. 
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1.3 THESIS LAYOUT 

In chapter 2, first an overview of some of the well established methods of stress 
separation is presented. It then discusses in detail the shear difference method and it’s 
modified version i.e., Trebuna’s method. Comparative analysis of the two methods is 
made by intentionally introducing error in the isoclinic parameter (0). The results have 
been graphically recorded and presented. The effect of grid size on the accuracy of the 
results has also been discussed. 

Chapter 3 focuses on the use of experimental data for stress separation. The 
analysis starts with theoretical solutions, being used as input data. Then partly realistic 
data is used for the analysis. This partly realistic data is in the form of the isochromatic 
parameter (N) which is obtained by the phase-shifting technique. Finally completely 
realistic data is employed for the analysis, with the isoclinic parameter ( 6 ), which in real 
sense is the main cause of errors, being experimentally determined using the circular 
polariscope arrangement and the plane polariscope arrangement. Here it is proved that 

0 obtained from plane polariscope arrangement is better than that obtained from circular 
polariscope arrangement. The superiority of 0 obtained by superimposition of two 

images (grabbed with plane polariscope arrangement) over 0 obtained by a single image 
is also discussed here. The analysis for the case of plate with a hole, using completely 
realistic data is also discussed in this chapter. 



CHAPTER 2 


STRESS SEPARATION TECHNIQUES IN 
PHOTOELASTICITY 


2.1 INTRODUCTION 


Photoelasticity is an experimental method for analyzing stress fields which yields the 
principal stress difference and its orientation in the form of fringes that are known as 
isochromatics and isoclinics respectively. 

The fringe pattern obtained by normal incidence in a circular polariscope can be 
related to principal stress difference using stress optic law as, 


a i ~°2 


NF^ 

t 


( 2 . 1 ) 


where F a is material stress fringe value (N/mm/fringe), N is fringe order at the point of 
interest and t is model thickness in mm. 

It is evident that the pure photoelastic method cannot solve completely a general 
stress problem. The conventional photoelastic method provides two independent 
equations for a two-dimensional problem and five independent equations for a three 
dimensional problem. Since there are three unknown stress components for a two 
dimensional problem and six unknown stress components for a three dimensional 
problem, the equations obtained by photoelasticity are not sufficient to determine the 
complete state of stress at any interior point. Hence with a view to obtain additional 
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relations, several auxiliary methods are required for both two and three dimensional 
problems, which in conjunction with photoelastic analysis can provide the complete 
information for solving a general stress problem. Determination of crl and cr2 separately 
by utilising additional information is referred as stress-separation. 

2.2 Methods of Stress-Separation 

In a broad sense, the stress separation techniques can be divided into four main 
categories. These categories are 

(i) Oblique Incidence Method. 

(ii) Use of Compatibility Equation. 

(iii) Shear difference Method or integration of the equilibrium condition. 

(iv) Hybrid methods based on combination of experimental & numerical methods. 

In oblique incidence method, additional equation for stress separation is obtained 
by making photoelastic measurements in oblique incidence, as shown in Fig. 2.1. To 
minimise refraction effects the model is to be immersed in a tank of liquid having same 
refractive index as the model. The secondary principal stresses along the light path are 
assumed to cause the formation of photoelastic fringes 4 . The method is basically a point 
by point technique. 

If the model is viewed normally( the direction of incidence being along z axis 
in x-z plane) and obliquely (the directions of incidence being at angles +(|) and -<|) to the 
z axis) the data obtained is sufficient to yield the unknown stress components. 
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Fig. 2.1 Oblique incidence method. 

Case a: Normal Incidence 

N : Corresponding Fringe order 

9 : Corresponding Isoclinic Parameter 
Case b: Oblique Incidence at angle <j>’ 

N’: Corresponding Fringe Order 

0’: Corresponding Isoclinic Parameter 
Case c: Oblique Incidence at angle -({)’ 

N": Corresponding Fringe Order 
9": Corresponding Isoclinic Parameter 

If the notation for the Fringe order and Isoclinic Parameter corresponding to the 
different directions of incidence, is as referred above and the material fringe value and 
the model thickness are represented by F c and h respectively, then the solution can be 


represented in the form: 
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N F 

cos26 

h 


( 2 . 2 ) 


F 

o -o = 2 — (N"cos2d" cos <p' + N'cos2d'cos(p' - 2 Ncos20) (2.3) 

2 hsin 2 <p 

t = — Nsin26 (2-4) 

^ 2 h 


FcOS(f> 

T = -2 (Nsin2d-N'sin2d') ( 2 - 5 ) 

K Iksmtp' 

F 

t = 2 — (N"cos26" -N'cos2&) ( 2 - 6 ) 

w 4sm0' 

This is the maximum information obtainable from photoelastic observations alone 
and enables the determination of the individual stress components. The biggest 
advantage of the method is that the stresses can be separated at all points in the model 
independent of neighbouring points. 

In the compatibility equation method, separation of stresses is accomplished by 
using the solution of ct x + o y , (the first stress invariant for constant body forces) and the 
photoelastic data a, - a 2 . This offers greater accuracy since the isoclinic angle 9 is not 
involved in the determination of principal stresses. The compatibility equation method 
requires a good knowledge of the boundary values of the model which may be difficult 
to achieve as time-edge effect or residual stresses due to machining are significant. 

Before the advent of widespread use of digital computers, the shear difference 
method was widely used technique for stress separation. One of the main drawbacks of 
this is the accumulation of error, which builds up as one proceeds along the line of 
interest. This has been reported mainly due to the poor experimental evaluation of the 
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isoclinic angle. The recent modification by Trebuna on numerically integrating the shear 
component’s differentials, has to a great extent diminished this problem. In conventional 
analysis the data is collected over discrete array of points. However, in digital 
photoelasticity, the data array could be finer and it can be just one pixel. The resolution 
of the pixel could be increased by appropriate optical magnification. The purpose of the 
present thesis is to have a relook on the applicability of shear difference technique for 
stress separation in digital photoelasticity. 

The hybrid techniques 5 are extremely involved. Some techniques utilize the 
photoelastic results partially and rely on numerical procedures while some attempt to 
obtain the required information based entirely on experimental data. The technique of 
using only the experimentally determined surface stresses and numerical methods for 
interior stresses facilitates important applications. These techniques have been utilized 
for stress analysis of 3-D bodies as well as on models of different material like photo- 
orthotropic materials 6 . 

2.3 NEED FOR APPLICATION OF SHEAR DIFFERENCE SCHEME IN 
DIGITAL PHOTOELASTICITY 

The main advantage of photoelasticity is that, it can give the stress distribution for the 
whole field. But it gives this result as an image of fringe patterns. The main task in 
photoelasticity is the interpretation of these fringe patterns. In conventional methods 
tedious compensation methods are needed to obtain the fractional fringe orders at any 
point other than the fringe areas. The number of data points that could be considered is 
also less. These procedures are quite involved and requires skill in the identification of 



10 


isochromatics. Hence, automation of data acquisition and analysis to minimize these 
problems and to provide fast and most accurate results, has become essential. With the 
advent of PC based digital image processing systems, considerable success has been 
achieved in automating the data acquisition from the entire field and it has now become 
possible to evaluate the total fringe order at every pixel in the domain. With the 
availability of photoelastic data, at pixel level, it becomes relevant to study shear 
difference scheme in the context of digital photoelasticity. This would offer scope for 
improvement in the results as unlike the conventional approach where a coarse data 
array is selected, the data array size could be much finer, as fine as just one pixel. 

2.4 SHEAR DIFFERENCE METHOD 


The shear difference method, uses the integral form of the equilibrium equation 7 . In the 
absence of body forces, the equation of equilibrium when applied to the problem of 
plane stress are, 


do* 

etc 

da 

y 

dy 


dr 


yx 


dy 

etc 


= 0 
= 0 


Equation (2.7) may be re-arranged and integrated to become. 




(2.7) 


( 2 . 8 ) 


(2.9) 


These can be approximated using the finite difference expressions, 
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= KX 


Ax 


( 2 . 10 ) 


°» - i5> 


( 2 . 11 ) 


In eq. (2.8) and eq. (2.9), (ctJo and (cr y )o represent known stresses at the 
boundary point chosen as the start of the integration procedure. If the starting point is 
chosen on a free boundary then (o x )o or (a y )o are given directly by the isochromatic 
parameter, is then calculated at the interior points and the summation proceeds in a 
stepwise manner using a fixed grid with each block measuring Ax by Ay. Ax xy is 
obtained from the difference between the top and bottom values at each interior point*. 
In the above equations, the convention for positive x is from left to right while for 
positive y is from bottom to top. Here in the present work, for the horizontal grid the 
direction of scanning is from outer (ie.boundary of the quadrant of the circular disc) to 
the interior (ie. line of symmetry). Since this direction of scanning ie. from right to left, 
is opposite to that of the standard convention as shown in Fig. 2.2, the sign of Ax needs 
to be monitored. With the sign of Ax being negative eq. (2.10) for o x takes the form 




( 2 . 12 ) 


Once a x is calculated for a particular point on the grid, a y can be calculated based 
on the relation : 


NF_ 


o = o 

y x 


cos20 


(2.13) 
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Fig. 2.2 Direction of scan along a horizontal grid. 

Hence cr x and <r y can be evaluated on each and every point of the grid, while 
scanning from right to left. It is to be noted that what shear-difference method demands 
is that the starting values be known. So, once the horizontal scanning is done and g x and 
a y at each point of the grid evaluated, these values can now be used as starting values 
for individual grids being formed for vertical scanning. While scanning vertically, first 
o y is evaluated and then is g x . The equations used are either of eq. (2.14) or eq. (2.15), 
depending on the location of the point of interest. 


- E 


At 


*y 


Ax 


Ay 


(2.14) 


+ E 


At. 


*y 


Ax 


Ay 


(2.15) 


Equation (2.14) is used when the points on the grid above the known point are 
being scanned while eq. (2.15) is used when those below the known point are scanned. 
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Once CT y is calculated ct x is calculated based on the relation 

NF n , 

o = a v + cos2f? (2-16) 

1 y t 

Likewise the stress values o x and a y are obtained for whole of the quadrant of 
the circular disc under diametral compression. 

2.5 TREBUNA’S METHOD 

This method suggests that the shear stress be expressed in terms of the difference of 
principal stresses, as 

Xfy = — Oj) sin20 (2.17) 

and differentiated it as a product of two functions. Substituting this into eq. (2.10) and 
eq. (2.1 1) yields 

a = (a) -f— — - — ——sm26dx-f(a.-a 1 )cos2d— dx (2-18) 

1 x o J dy 2 J 1 2 dy 

a y = ( a A"/“^aT^i sil120 dy-f(o r o 2 )cos2d^- dy ( 21 9) 

These can be approximated using the finite difference expressions as, 

°x = KXrE - A( — — | sin2g A *-£ (°r^ cos2 ^ A * (220) 

° y = («Po-E Ay-^(a r a 2 )cos20-^ A y ( 2 -21) 


Again, as for the horizontal grid the direction of scanning is opposite of the 
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standard sign convention, the sign of Ax needs to be monitored. Hence a x is now 
approximated as 


°jc = Ax+^(aj-a 2 )cos20-^ Ax 


( 2 . 22 ) 


= 


NF 

cos20 

t 


(2.23) 


Hence a x and 0 y can be evaluated on each and every point of the grid, while 
scanning from right to left. Now again these points act as starting points for vertical 
scanning.While scanning vertically we first evaluate a y and then a x .The equation used 
is either of eq. (2.24) or eq. (2.25), depending on the location of the point of interest. 


* K).-E —.LtJ istize Ay-£ (o,-a 2 )cos2e4^ Ay (224) 

y r Ax 2 Ax 

°y = ~T — — ^-sin2(9 Ay+J^ (o r o 2 )cos20— Ay ( 2 - 25 ) 

y r Ax 2 Ax 

Equation (2.24) is used when the points on the grid above the known point are 

being scanned while eq. (2.25) is used when those below the known point are scanned. 

Once cj y is calculated, cr x is calculated based on the relation 



+ 


NF„ 

cos20 

t 


(2.26) 


2.6 GENERATION OF THEORETICAL DATA AND APPLICATION OF 
STRESS SEPARATION SCHEMES 


The earlier section discussed the standard form of equations related to shear-difference 
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method & Trebuna’s modification to the shear difference technique. It also discussed 
their modified form, in case the direction of scanning is opposite to the conventional +x 
or +y axis. This section discusses the implementation of these equations over the 
geometry under consideration i.e., disc under diametral compression. In the present 
analysis only a quarter section of the disc is considered and the whole field analysis of 
the quadrant is desired. 

2.6.1 Physical significance of the grid 

One of the essential feature of both the above mentioned methods is that they need to 
be provided with the starting values of a x or a y referred to as (a x ) 0 Sc (cr y ) 0 respectively 
so that the integration procedure can be applied. 

In the case of disc under diametral compression; the isochromatic parameter ie 
the fringe order on the free boundary is zero & hence ( a x ) 0 or (cr y ) 0 can be directly 
obtained. With the starting point taken on the boundary, we can scan from boundary to 
the line of symmetry and obtain the state of stress at all points located at specified 
regular intervals. This interval is at the discretion of the user. The smaller the interval, 
the better the integration procedure & hence better will be the results. It is noteworthy 
that these points with known state of stress can now act as the initial point for vertical 
scanning i.e., each of these points will act as the starting point for the integral procedure 
in both +y Sc -y directions, so as to obtain the results over the whole-field. 

It’s also important to note from eq. (2.8) and eq. (2.9) that for shear-difference 
method, the variation of shear-stress is averaged over a specified increment along x or 
y direction. Similarly, in Trebuna’s method the variation of (cr ; - <s 2 ) (which implies the 
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variation of isochromatic parameter N) and the variation of isoclinic parameter 0 is 
averaged over a specified increment along ror; directions, as reflected by eq. (2.18) 
& eq. (2.19). This physically emphasises the need for grids, wherein the centre line is 
the line of interest & the adjacent ones are to facilitate the variation of the desired 
parameter along x or y directions. Figure 2.3 shows the grid system, to be adopted. For 
better understanding, consider eq. (2.8) for the case of shear-difference method, a x at 
point B would be given by 

K)a = (°X ~ delx (2.27) 

dely 



Fig. 2.3 Concept of a grid system. 
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Similarly, consider eq. (2. 18) for the case of Trebuna’s method, wherein cr x atthe 
point B would be given by 


(° j)b = (°*>o - 


( a r a 2 >A - (°f° 2 )c 1 


(2.28) 


(grgj 

dely 


defy 


sin2d D delx 


(Oj -o 2 ) s cos26 g delx 


The same concept of grid, is applicable for vertical scanning also. 


2.6.2 User’s input data 


In order to obtain the whole field state of stress for the quadrant of the disc under 
consideration, the pre-requirement is 

(i) a known, diameter of the disc. 

(ii) a known, line of interest (yd). 

(iii) a known, increment along x-axis. 

(iv) a known, increment along _y-axis. 

These parameters in fact constitute the data which is at the discretion of the 
analyst. Under-given is a detailed procedure with regard to the shear-difference method 
of stress-separation. It remains exactly the same for Trebuna’s method except for the 
applicable equations. 


2.6.3 Scheme for horizontal scanning 
(i) No. of data points: 

Once the above mentioned set of data has been fed in, the leng ' line of interest 
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along x axis say xd can be known by the relation : 

xd = sqrt(/^ 2 -(yd) : ) 

Then the no. of divisions along the line of interest = xd/delx. Based on this value the 
no. of data points along the line of interest say ndp are evaluated. Two separate cases 
thus arise, which are shown in Fig. 2.4. 
case a) All grids of equal size. 

eg. xd = 10.4 mm; delx = 0.52 mm 
no. of divisions = 10.4/0.52 = 20 
ndp = 21 ( each at equal interval ) 

This is the case wherein all the grids will be of equal size. 



Fig. 2.4 


Scheme for horizontal scanning: (a) All grids of equal size (b) Unequal 
grid at boundary. 
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case b) Grids of non equal sizes. 

eg. xd = 10.7 mm; delx = 0.52 mm 
no. of divisions = 10.7/0.52 = 20.576 

ndp = 22 (21 at equal interval; one at non-equal interval = 0.576 times 
of delx ) 

diffx = 0.576 * delx. 

This is the case wherein all grids except one, will be of equal size. 

(ii) Generation of co-ordinates: 

For any photoelastic analysis over a domain, the x & y co-ordinates, the 
isochromatic parameter N and the isoclinic parameter 9 are necessarily required over 
each point of interest. The generation of x, y co-ordinates is discussed here, 
case a) The coordinates are to be generated over 3 lines. 

(i) A line, over the line of interest, referred to as Top line 

(ii) line of interest 

(iii) A line, below the line of interest, referred to as Bottom line. 

As the co-ordinate of the point on the free boundary on the line of interest is 
known from the relation xd = sqrt(R 2 -(yd) 2 ), the 0 th point of the top line can be 
allocated as it’s x coordinate the same value as xd and then each time delx can be 
deducted so as to yield x-coordinates till line of symmetry. 

It’s true that 0 th point of the top line, is rather in free space and not exactly on 
the free boundary but this is tolerable pertaining to the fact that the top line is just a 
pixel increment away from the line of interest and hence the error would just be 
negligible. The y-coordinates for all the points of the top line = yd + dely. Similarly the 
x-coordinates for the line of interest & bottom line can be allocated. For these two lines, 
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y-coordinates differ, as , for all the points on the line of interest y = yd; while for 
bottom line y = yd - dely. 

case b) In this case, the non-equal increment referred as diffx is to be accounted. It’s 
a better choice to account it near the boundary owing to the fact that the integration 
procedure starts from the boundary and hence integration performed over a small 
increment at an initial stage would yield a better solution. The generation procedure for 
all points remains exactly the same as case a except for the points 1, ndp+1, 2ndp+l 
whose coordinates equals to x boundary - diffx and not x bouildiry -delx. 

(iii) Integration Procedure. 

case a) a x for all the points from l-j can be determined from eq. (2.12) while cr y from 
eq. (2.36). 

case b) Determination of ct x is slightly different for point no. 1 For this Ax of eq. (2.12) 
is to be replaced by diffx. Once a x is determined, a y can still be calculated from eq. 
(2.36). 

With the procedure discussed above, being followed, the results obtained for 
three different scan lines, each measuring approximately y = 0, R/2, 3R/4, respectively 
have been presented in Fig. 2.5(a), 2.5(b) and 2.5(c) respectively. The results show that 
the stresses obtained by shear difference scheme and Trebuna’s scheme are in close 
agreement. In fact the stresses evaluated from the two methods, for the line y = 0 (mm), 
exactly overlap with the stresses evaluated, using the standard formulations for the disc, 
based on theory of elasticity (refer appendix-I). However slight variation can be 
observed for the line y = R/2, which further aggravates slightly for the line y = 3R/4. 
Even for these two lines, the results obtained by shear difference and Trebuna’s scheme 
overlap, though they vary slightly with the theory’s. 
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2.6.4 Scheme for vertical scanning 

(i) No. of data points: 

By this stage, a x and a y at all desired points along the horizontal line of interest are 
known. Consider a general point i and say it’s x coordinate is vxd. Now the length of 
the vertical line of interest can be determined by the relation vyd = sqrt(R 2 -(vxd) 2 ). 
Two separate cases again arise. Fig. 2.6 physically explains these cases. 



Fig. 2.6 Scheme for vertical scanning: (a) All grids of equal size (b) Unequal grid 
at boundary. 
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case a) All grids of equal size. 

eg. vyd = 19 mm 

dely = 0.76 mm 

no. of divisions = 19/0.76 = 25 

ndp = 26 ( each at equal interval ) 

This is the case wherein all grids will be of equal size, 
case b) Grids of non equal sizes. 

eg. vyd = 19.5 mm 
dely = 0.76 mm 

no. of divisions = 19.5/0.76 = 25.657 

ndp = 27 (26 at equal interval; one at non-equal interval=0.657 times of 
dely ) 

diffy = 0.657 * dely 

This is the case wherein all grids except one, will be of equal size. 

(ii) Generation of coordinates: 

The procedure of generation remains exactly the same as for the horizontal 
scanning. The lines are referred as Right line, line of interest & the Left line. The x 
coordinates for all points along : 

Right line = vxd + delx 
Line of interest = vxd 
Left line = vxd - delx 

The y coordinates are allocated from the boundary towards the line of symmetry (x- 
axis). For case(a) mentioned above, this allocation is uniform while for case(b). diffy 
is to be accounted close to the boundary. 
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(iii) Integration procedure: 

The integration procedure for vertical scanning is relatively more involved. Four 
considerations are to be made. First two corresponding to the cases (a) and (b) as 
mentioned above, i.e, the grids being of equal or non equal size; while the next two 
corresponding to the assignment of the values of the known points on the horizontal line 
of interest as the starting values for vertical integration for each individual grid. 

It’s worthy to realize that the dist. yd fed in by the analyst may not necessarily 
be an integral multiple of dely. So, the direct assignment of the known values of the 
horizontal points as the starting values for vertical integration may not be possible; 
forcing two new aspects into consideration. 

Aspect(i) 

yd = n * dely ; where n : integer 

As shown in Fig. 2.7(a), the known point say ’k’ (known from the horizontal scanning) 
by virtue of it’s overlap as one of the point of interest along the vertical line; can be 
directly used as the starting point for the integration procedure. Now with the starting 
value known, utmost care should be exercised while determining the applicable 
equations for calculation of a x & cr y for points below and points above the known point. 
As mentioned and justified in earlier section eq. (2.14) & eq. (2.16) are valid for points 
above while eq. (2.15) & eq. (2.16) for points below the known point. 

Note: a 0 = cr k 

Aspect(ii) 

yd != n * dely ; where n : integer 

Fig. 2.7(b) shows that the known point ’k’ (known from the horizontal scanning ) does 
not overlap with the point of interest along the vertical line, and hence it cannot be 
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directly used as the starting point for the integration procedure. Rather, the known value 
is to be transformed to the point of interest so that it may act as the starting point. This 
transformation is not done arbitrarily or by approximation but by applying the equations 
proposed in the shear difference & Trebuna’s method for vertical directions. As a 
predetermined convention, the known values are transformed to a point of interest which 
is just below the known point & also happens to be at a dist. which is an integral 
multiple of dely. 

For better illustration of this transformation, consider the known point k and the point 
of interest which also has to be the starting point. The transformation is as follows: 

(op, - (op 4 * (T>A (fe k (T ' j) ‘ de ‘y (2 29) 

a x can be calculated using eq. (2.16). 

Note: a 0 = a. 

Now again as the starting values are known, attention must be focused on the 
permissible equations corresponding to : 

- different scan directions ( above or below known point ) 

- two different cases: case a) and case b) respectively, 
case a): 

For evaluation of points below the starting point, eq. (2.15) & eq. (2.16) are 
valid, while eq. (2.14) & eq. (2.16) are applicable for points above starting point, 
case b): 

As for case a) eq. (2.15) & eq. (2.16) provide the stresses for points below the 
starting point while for points above the starting point, the strategy is slightly different. 
Equation (2.14) & eq. (2.16) are applicable till the point just below the free boundary. 
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but for calculation of the point on the free boundary, the term dy in eq. (2.14) for 
relation of cr y is to be replaced by diffy; diffy being the non equal increment accounted 
at the boundary. 

The discussion over the vertical scanning, conveys the sense that once the 
horizontal scanning is done and stresses at the points along this line obtained, then each 
of these points can act as the initial point for vertical scanning. Depending on the 
location of the horizontal line of interest, the scanning in vertical direction can be 
performed both along the +y and -y directions to cover the region of interest. Fig. 
2.8(a), 2.8(b) and 2.8(c) present the results for three vertical lines, measuring 
approximately x = 0, R/2, 3R/4 (mm) respectively, with the horizontal line at y = R/2, 
taken as reference. The results for the horizontal line at y = R/2 have already been 
presented in Fig. 2.5(b). Along this line, the points at x = 0, R/2 and 3R/4, with known 
stress values were taken and integrated in +y and -y directions according to the shear 
difference scheme and the Trebuna’s scheme. It can be observed that the results 
produced by these methods are ic close aggreement with the theory’s. 

The complete procedure mentioned so far was for a general point, known over 
the horizontal line of interest (implying a general case of a vertical grid). The same 
procedure can be applied to all the known points over the horizontal line of interest, thus 
enabling the whole field analysis of the disc under diametral compression. Figure 2.9 
shows the field variation of the stresses. Figure 2.9(a) and Fig. 2.9(b) represent the 
variation of ct x and a y respectively, evaluated by shear difference scheme. It has been 
concluded from the results for various horizontal scan lines that with the theoretical data 
being employed, the results obtained by the standard solution for the disc under 
diametral compression, the shear difference method and Trebuna’s method for stress 
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separation, are in complete agreement. As, the field variation of stresses obtained by 
theory and by Trebuna’s method are exactly identical to those of shear difference 
method, they have not been shown explicitly. 

2.7 EFFECT OF GRID SIZE 

In this present approach, which aims at the establishment of shear difference and 
Trebuna’s method as reliable stress separation schemes for whole field analysis of 
stresses, the required photoelastic data is theoretically generated by using standard 
formulations based on theory of elasticity (refer appendix I). Later, the use of 
experimentally obtained photoelastic data, is discussed. For this purpose, during 
experimentation, the images of the model under analysis are captured and the intensities 
of the recorded images are processed. The fact that the intensity information is processed 
for each pixel in the region of interest, limits the grid size to the pixel increment only. 
However there is no such restriction in this approach as the photoelastic data can be 
theoretically generated for any grid size. Fig. 2.10 and Fig. 2.11 represent the effect of 
grid size over the accuracy of the seperated stresses. Fig. 2.10 presents the results for 
horizontal scan line at y = R/2. It can be observed from Fig. 2. 10(a), 2. 10(b) and 2. 10(c) 
that the accuracy of the separated stresses goes on reducing as the grid size is 
sequentially raised from pixel dimensions to five times the pixel dimensions and finally 
to ten times the pixel dimensions. This observation is further supported through the plots 
in Fig. 2.11 which present the results with varying grid sizes for vertical scanning along 
the line at x = R/2 with the horizontal line at y = R/2, taken as the reference. It can be 
concluded that, the smaller the grid size, the better are the results. 









Magnitude of strata (MPa) Magnitude of strata (MPa) Magnitude of stress (MPa) 
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Effect of grid size over accuracy of stresses for vertical scan at x = R/2 
(mm), for an initial horizontal scan at y = R/2 (mm), (a) Pixel 
dimensions (b) 5 times pixel dimensions (c) 10 times pixel dimensions. 
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Fig. 2.12 Comparitive study of Trebuna vs Shear difference method for horizontal 
scan at y = R/2 with erroneous Isoclinic data: (a) 10% error (b) 20% 
error (c) 50% error. 







Fig. 2.13 Comparitive study of Trebuna vs Shear difference method for vertical 
scans with: 20% error in isoclinic data. Initial horizontal scan at y= R/2. 
(a) x = 0 (mm) (b) x = R/2 (mm) (c) x = 3R/4 (mm). 
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2.8 COMPARATIVE ANALYSIS OF ACCURACY ASPECTS OF SHEAR 
DIFFERENCE AND TREBUNA’S METHOD 

Trebuna’s method is supposed to be more accurate as the isoclinic parameter is averaged 
along the three lines of the grid over the model and assumed to represent the isoclinic 
along the centre line. The purpose is to "smooth out" any noisy points of isoclinic data. 
Here, in the present approach, as the isoclinic data is theoretically evaluated it is already 
good and hence its averaging does not introduce any improvement in (Trebuna’s) the 
solution over the shear difference solution. When error in isoclinic parameter is 
intentionally introduced and trebuna’s method then employed, it yields a solution which 
is more accurate than the one attained by shear difference. This effect has been reflected 
by Fig. 2.12 and Fig. 2.13. It is observed that with the isoclinic parameter turning 
erroneous, the results obtained by Trebuna’s scheme are more accurate than those 
obtained by shear difference scheme. Fig. 2. 12(a), 2. 12(b) and 2.12(c) reveal that, as the 
error in the isoclinic parameter goes on increasing, the level of additional accuracy that 
trebuna’s modification offers, also goes on increasing. Fig. 2.13 in support of the above 
observation shows the results for three different vertical scan lines, for the moderate 
value of error in the isoclinic parameter. 



CHAPTER 3 


STRESS SEPARATION IN DIGITAL 
PHOTOELASTICITY 


3.1 INTRODUCTION 

The recent advances in digital photoelasticity 9 have offered promise for full automation 
of data acquisition in two dimensional photoelasticity. Following this, a new whole field 
technique for determination of isochromatic and isoclinic parameters using phase shifting 
methodology is discussed and implemented. In this technique, a few images 
corresponding to different optical arrangements, are to be captured. Each of these 
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different optical arrangements, introduces a specific phase- shift and hence the technique 
is referred as phase shifting. Phase-shifting generally provides the value of isochromatic 
parameter (TV) with a high degree of accuracy. However,' the determination of isoclinic 
parameter (0) by the existing methods needs to be further improved for its practical 
utilisation. To determine (Q), two different optical arrangements are discussed. One 
involves the use of a circular polariscope arrangement and the other, the plane 
polariscope arrangement. 

In this chapter, the use of shear difference technique and Trebuna’s method for 
stress separation is discussed in thecontext of digital photoelasticity. Initially, the stress 
separation algorithm is verified by using the theoretical values of (N) and (6) for the 
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problem of a disc under diametral compression. Later stress separation is effected using 
the data of N obtained by phase shifting and (6) by theory. Finally the usefulness of the 
technique by using the data fully from experiment is investigated. 

3.2 NEED FOR AUTOMATION OF DATA ACQUISITION 

The main advantage of photoelasticity is that, it can give the stress distribution for the 
whole field. But it gives this result as an image of fringe patterns. The main task in 
photoelasticity is the interpretation of these fringe patterns. In conventional methods, 
tedious compensation techniques are needed to obtain the fractional fringe orders at any 
point other than the fringe areas. The number of data points that could be considered is 
also less. These procedures are quite involved and require skill in the identification of 
isochromatics. Hence, automation of data acquisition and analysis to minimize these 
problems and to provide fast and more accurate results, has become essential. With the 
advent of PC based digital image processing systems, considerable success has been 
achieved in automating the data acquisition for the entire field. In the early development 
of digital photoelasticity, the data acquisition is limited to fringe areas. The real potential 
of DIP hardware was realized only when the concept of identifying the fringe fields as 
phase maps came into existence. This introduced new concepts in data acquisition. It is 
possible to evaluate the total fringe order at every pixel in the domain and not restricted 
to just fringe skeletons alone. 
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3.3 Phase shifting In photoelasticity 

The phase stepping concept can be simply described as "changing the absolute phase of 
the reference wave in equal steps and measuring the local light intensity after each step". 
The change in phase, is in practice, achieved by rotation of the output optical elements 
of the polariscope. By varying the phase difference between the beams involved, in 
known steps, it is possible to generate sufficient number of equations to solve for the 
parameters involved. The phase shifting approach has been applied to many classical 
interferometric techniques. In most classical interferometers including holography, the 
path length of the two interfering beams are distinct and separate. In general, phase 
differences can be added by altering the path length of any one of the light beams. In 
photoelasticity, the two light beams cannot be treated separately but always go together. 
This means that a phase shift introduced in one light beam will also introduce a 
corresponding phase shift in the other beam also 10 . The change in phase, in practice, is 
achieved by appropriately rotating the optical elements of the polariscope. 

The phase shifting concept for photoelasticity was first introduced by Hecker and 
Morche 11 in 1986. Extending the work of Hecker and Morche, Patterson and Wang 12 
reported in 1991 a full automation procedure for photoelastic analysis. In 1992, Sarma 
et al. 13 reported a methodology using plane polariscope. Later in 1993, Asundi 10 reported 
a new methodology, which extended the Tardy method of compensation from a point 
to point approach to evaluate fractional fringe order on all points lying on the isoclinic 
line. In 1996, Ramesh and Ganapathy 14 reported a comparative study of the above 
methodologies and reported that among these techniques, the technique proposed by 
Patterson and Wang is truly a whole field technique. 
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Fig. 3.1 Circular polariscope arrangement for phase shifting technique. 
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3.3.1 Patterson and Wang’s Scheme;Use of circular polariscope. 


The Patterson and Wang’s algorithm requires a circular polariscope which has a 
provision to rotate each optical element independently. When circularly polarised light 
is used to view a photoelastic specimen, the light emerging from the analyser at any 
point (x,y) takes the form for intensity: 

i(x,y). = i m + i v (sin2(/3-<p) cos<5 - sin2(d-tp ) cos2(fi-(p) sin<5) (3.1) 

where i v is the intensity emerging when all the axes of the polariscope and specimen are 
parallel and the term i m takes account of the stray light. The angles /3, cj> and 6 are angles 
between a reference axis and the slow axis of the analyser, output quarter wave plate 
and specimen respectively, as shown in Fig. 3.1. The relative phase 5 of the two rays 
formed when the light is transmitted through the specimen is related to the fringe order 
N as follows: 


_6_ = N _ (°i ~ fa) t 

2n ~ ~ / ' 


(3-2) 


where a, and a 2 are the principal stresses in the plane of the specimen, t is the thickness, 
and / the material fringe constant. In eq. (3.1) there are usually four unknowns, namely 
the isoclinic parameter (8), the relative phase 8 and the intensities i m and i v . In order to 
solve for these unknowns a minimum of four phase steps or changes in 0 and <f>, the 
positions of output optical elements, are needed. Hecker and Morche employed 
increments of ir/3 for /3, the position of the analyser with the output quarter plate both 
aligned to the analyser and at ± x/4 to it. They suggested that the use of high 
modulation of the phase produces smoother results. Consequently Patterson & Wang 
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used increment of 7 t/ 4 for both (3 and <£, and combined the resulting six signals in such 
a way that a modulation of x/2 is achieved. 

Below, are presented in a tabular form, the various possible combinations of 
Analyser angle (3 and the angle <j>, representing the orientation of the second quarter 
wave plate. The intensity equations corresponding to each combination of 13 and 4> are 
also specified. 


Table 3.1 Intensity equations for various optical arrangements used in the phase- 
shifting algorithm. 


S.No 

II Quarter 
wave plate 

Analyser 

angle 

P 


Intensity equation 

Fringe pattern 
interpretation 

1 

0 

7t/4 

h 

= i' + i L cos 6 

tn v 

circular polariscope 
bright field 

2 

•0 

3tt/4 

h 

= i - i„ cos 6 

m v 

circular polariscope 
dark field 

3 

0 

0 

h 

= / - i„sin2 0sin 8 

in v 

none 

4 

7l/4 

7l/4 

*4 

= i‘ + i cos2 6>sin(5 

m v 

none 

5 

tc/2 

n/2 

h 

= L + ism26sm6 

m v 

none 

6 

3tc/4 

3tc/4 

h 

= i - i cos20sind 

m v 

none 


Using the above intensity equations, the isoclinic angle ( 8 ) and the fractional retardation 
(6) can be evaluated as mentioned below. 
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g. 3.2 


Plane polariscope arrangement for phase shifting technique. 
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3.3.2 PLANE POLARISCOPE ARRANGEMENT 

In phase shifting technique, the relative position of the optical elements play a 
significant role. Removing the first and the second quarter wave plate from the circular 
polariscope arrangement and changing the relative orientation of the analyser and the 
polariser 15 , leads to a new arrangement. As only the polariser, the model and the 
analyser are involved in the arrangement, it is referred as plane polariscope arrangement 
(refer Fig. 3.2). If the orientation of the polariser is denoted by angle (a) while that of 
the analyser by angle (/3) and the intensities corresponding to these combinations 
recorded, it is observed that only six intensity equations turn out to be unique for all 
possible combinations of a and /3, in steps of 7r/4. 

Table 3.2 Table, showing that only six intensity equations are independent. 
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These six intensity equations corresponding to various combinations of the a and /3 have 
also been formulated in table 3.3. Using the above mentioned intensity equations, the 
isoclinic parameter (0) can be evaluated using the relation: 
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0 = — tan 


1 . h~h 


-i 2 . 2 

i *4 " *1 


Table 3.3 Table showing the intensity equations for the different orientations of the 
polariser and analyser. 
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Polariser Analyser 
angle a angle p 
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The fractional retardation can still be evaluated using the eq. (3.4). 

3.4 Image Processing Set-up 

The image processing system consists of a CCD (charge coupled device) camera (TM 
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560 PULNIX) that has a pixel resolution of 512 x 512 pixels and digitizes the images 
at video rates. 1 he camera is connected to a PC based image processing system 
equipped with an image processing card (PIP 1024B, Matrox corporation) connected to 
a high resolution video monitor. The system has in addition to PC memory, four frame 
buffers each of 512 x 512 pixels x 8 bits, to store the image being processed and the 
intermediate results. 

3.5 Experimental Evaluation and Whole Field Representation of 

Fractional Retardation 

For the above mentioned experimentation process Sodium vapour lamp is used as 
monochromatic light source. An epoxy disk with a diameter of 60.06 mm, thickness of 
5.04 mm and material stress fringe value (F g ) of 1 1.64 N/mm/fringe, which is subjected 
to a diametral load of 502.57 N is taken up for the analysis. Initially, adopting tire 
circular polariscope arrangement, six images are recorded by keeping the second quarter 
wave plate and the analyser in six different positions as shown in Table 3.1. Similarly, 
images are also captured for six different optical positions as proposed in plane 
polariscope arrangement. These optical positions are referred in Table 3.2. The captured 
images are shown in Fig. 3.3. The eq. (3.3) produces values of (8) which are periodic 
with a wave length x/2. 

It is obtained in a range from -7r/4 through zero to 7r/4 which corresponds to the 
range of 0 to 7r/2, obtained from the classical definition of the isoclinic angle given by 
Frocht7. The value of the phase retardation is expressed in the range of -x to x by using 
atan2 function of C-program library. Dark field isochromatics can be obtained using the 
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Fig. 3.3 The sequence of six images of a circular disc under diametral 
compression recorded for phase shifting technique, by plane polariscope 
arrangement. 
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following equation. 

g (x,y) = abs( - 6) - 255.0 (3-5) 

TC 

Where g(x,>') is the grey level of the pixel at the location 



Fig. 3.4 Figure showing the various representations of the phase retardation. 

(a) Simulated dark field (b) Simulated bright field (c) Phase map. 

Figure 3.4(a) shows the simulated dark field isochromatics. The phase map is obtained 
by converting phase retardation to grey levels. The zero phase retardation corresponds 
to pitch black and 27r phase retardation corresponds to pure white. So, the phase 
retardation is to be expressed in the range of 0 to 2 tt. The phase map can be obtained 
from the following equations. 




5 = 5 
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5 = 5 + 2r 8 < 0 

g(x,y ) = 255.0 (3.7) 

In 

Figure 3.4(c) shows the phase map of the image. The bright field isochromatics for the 
image can be obtained using the following equations. 

5 = 5 + 2tt 5 < 0 

g(x,y) = 255.0 (3.8) 

71 

Figure 3.4(b) shows the simulated bright field isochromatics of circular disk under 
diametral compression. 

3.6 Boundary Detection 

To limit the determination of field parameters within the region of interest (within the 
specimen), the image has to be segmented. The segmentation of the image will result 
in isolation of the region of interest from the rest of the image. Segmentation of the 
image is done by detecting the boundary. In general, boundary of the image is detected 
using spatial filters like Roberts, Prewitt, Sobel etc.. These filters detect the edges based 
on abrupt change in grey levels. The edge detection algorithms not only detects the 
boundary, but also detects the regions in the image where there is abrupt change in grey 
levels. Moreover the primary interest is the detection of the boundary and not the 
detection of edges. The set of pixels obtained from above filters seldom characterizes 
a boundary completely, because of noise and breaks in the boundary from non uniform 
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illumination. These edge detection algorithms typically are followed by other boundary 
detection procedures designed to assemble and link the edge pixels into meaningful 
boundaries. This only yield raw data in the form of pixels along a boundary. These 
filters are helpful to detect the edges of the images obtained in DIP applications like 
astronomy and defense where image recognition is the only criterion. Accuracy of the 
boundary is not critical. Unlike the above applications, the experimental mechanics 
application like photoelasticity, the shape and boundaries of the object are known. This 
property is used to detect the boundary accurately instead of using the above general 
edge detection algorithms. To make computations, the data is to be expressed as a 
suitable representation. In the present study the algorithm proposed by Sreedhar is used. 
In this the boundary is identified by defining two files one for row-scan and the other 
for column scan. The area of the analysis is limited by drawing the boundary around the 
region of interest. In the bright field arrangement of polariscope, the boundary of the 
image can be seen clearly [Fig. 3.3(a)]. To mark the boundary, primitives have been 
developed Sreedhar 16 . These are line, circle, ellipse and 3-point arc. The boundary is 
drawn with an intensity of 255. 

In order to quicken the boundary detection process, the algorithm accepts the 
actual dimensions of the specimen if available. For example to identify the circular disk 
boundary, one has to specify the centre and the radius. Based on the optical 
magnification and aspect ratio of the video monitor the boundary is drawn. Figure 3.5 
shows the histogram of the bright field image of the circular disk under diametral 
compression. It can be observed that the maximum grey level of the image over the 
entire image is 168, which is less than 255 grey level. So, any pixel with grey level 
value of 255 represents the boundary pixel over the entire domain. 
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Histogram of the bright field image of the circular disc under diametral 
compression. 



Magnified pixel representation of a curved segment on the screen 
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Boundary detection for a circular disc 





The image is represented as a grid of pixels and the representation of a curved 
boundary is shown in Fig. 3.6. In view of spatial quantization, to represent tire curvature 
more than one pixel is needed. However, for all computational purposes only the 
outermost pixel must be recognized as a boundary. In order to simulate the actual 
boundary, the pixel outside the actual boundary but just enclosing it should be marked 
as the boundary pixel. The boundary is drawn using the primitives developed such that 
the outer side pixels of the circle represents the boundary of the circular disk. In the 
ideal case, for the curve shown in Fig. 3.6 only one pixel should represent the boundary. 
This is achieved by storing the positions of the pixels which have an intensity of 255 
into two separate intermediate files (scanx.dat and scany.dat), one for horizontal 
scanning and the other for vertical scanning. These files will have more than one pixel 
forming the boundary. Further processing of these files is done in which the first and 
last columns for every line from these files are extracted and stored into two separate 
files for horizontal scanning ( *.ybn ) and vertical scanning (*:xbri); 

In the present study, the analysis of the image is limited to first quarter (256 x 
256) of the quadrant. Since, the specimen will not occupy the full quarter of the 
quadrant, some blank lines appear in scanned files. The boundary pixels corresponding 
to these blank lines are 0 and 255. Hence, these blank lines are replaced with 0 and 255. 
The boundary detected by applying the above method on the disk is shown in Fig. 3.7. 


3.6.1 UTILITY OF BOUNDARY INFORMATION 


central umm * 

1 1 t., 

A 128076 


It has already been mentioned thus far, that two separate files, one each for horizontal 
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scanning (* .xbri) and vertical scanning ( *.ybn ) are formed to store total boundary 
information. Consider a disc being scanned over a window size of 256*256 pixels, with 
the left top comer of the window taken as reference (0,0). Consider two separate lines, 
line-1 and line-2 as shown if Fig. 3.8(a). Let line-1 intersect the region of interest at 
points A and B characterized by points P, and P 2 respectively. P, and P 2 represent the 
no. of pixels, that points A and B are away from the reference comer (left-top) of the 
window, along horizontal direction. Similarly, corresponding to line-2 let the points A ' 
and B’ be characterized by Q, and Q 2 , which denote no. of pixels along vertical 
direction. The values P, and P 2 are known from ( *.ybn ) file while Q, and Q 2 from 
(*.xbri) file. The aim lies in determining both the x and y co-ordinates for any general 
point along the boundary, which is possible through combination of information 
obtained by (*.xbn) and (*.ybn) files. In this regard if points A and A ’ are made to 
coincide, (Fig.3.8(b)) then the point of intersection would be characterized by (P,, Qi). 
Similarly let the centre of the disc be characterized by ( r,s ). Utilizing the symmetry of 
the disc, only a quarter section can be analyzed. For ease of analysis the centre of the 
disc must be regarded as the new origin (0,0). Then the intersection point A-A ’ will be 
characterized by (Pj-r^-Q,). The actual location of the point A-A ’,in terms of spatial co- 
ordinates can be given by (x,.y) where 

x (mm) = (P,-r) * delx 
y (mm) = (s-Q,) * dely 

Figure 3.9 represents the digitized boundary for the quarter of the disc i.e., the 
boundary determined based on the spatial location of the pixel on the boundary. It is to 
be noted that the (*.xbn) and (*.ybn) files determined the boundary in terms of the 
number of pixels that it is away from the reference taken on the video monitor, which 
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is the left top corner of the window containing the image. With the transformation of 
this reference to the centre of the disc and combining of the information given by 
(*.xbn) and (*.ybn) files, as mentioned above, the boundary for only the quarter of the 
disc has been determined and presented in Fig. 3.9. Similarly for all points along the 
boundary, the spatial co-ordinates are achieved. These known co-ordinates along the 
boundary help in determination of the co-ordinates of any general pixel inside the region 
of interest. Thus the boundary information enables the determination of spatial co- 
ordinates of all points in the region of interest. 

3.7 The grid system in digital domain. 

With photoelastic data available for each pixel in the region of interest, from different 
sources as previously discussed, stress separation is performed by scanning the image 
row by row from top to bottom. The algorithm developed for the implementation of 
shear difference and Trebuna’s Method, visualises the model as being constituted of 
three separate parts; the top-most part, the central part and the bottom-most part. The 
implementation of shear difference method and Trebuna’s method requires formation of 
grid constituted of three adjacent rows of pixels, with the central row being the row of 
interest i.e., the row corresponding to each pixel of which, the stress in x and y direction 
will be predicted. For the top-most and the bottom-most row of pixel, formation of an 
exactly identical grid may not be possible due to the absence of physical geometry 
above the top-most row and below the bottom-most row. Hence to perform the stress 
analysis over these rows, the averaging scheme of stress separation methods under 
consideration needs to be modified. Apart from the missing rows above or below, even 
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with adjacent rows available on both sides difficulty may arise due to non equal number 
of pixels amongst the available rows themselves. For better understanding consider a 
fully populated Grid-system as shown below. Referring to the eq. 2.10, eq. 2.11, eq. 
2.20 and eq. 2.21, it can be realised that with stresses known at the starting point and 
photoelastic data available for all other points, the determination of stresses for other 
points along the direction of scan requires the variation of shear stress or the variation 
of principal stress difference, to be known in a direction perpendicular to the scan 
direction. For ease of analysis, refer Fig. 3.10. 



Fig. 3.10 


Averaging scheme for a particular grid. 





(a) TOP MOST GRID 



(i) a 


(ii) b 


(iii) c 


(c) BOTTOM MOST GRID 
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Fig. 3.11 The grid system in digital domain. 
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3.11(a). Hence the parameter (t x) b~( 7 x) c f° r shear difference method while the 
parameters (o i~a b -(o i~a c and (0) b -(6) c for Trebuna’s method are required. Now, these 
need to be averaged only over one pixel increment along Y axis. 

CENTRAL GRIDS : 

Grid structure: All Rows (i), (ii) and (iii) available. 

Row of interest: Row (ii) 

Averaging Rows/Dist.: (a) For points (1) to (i) : Averaging rows (ii) and 

(iii); Dist=dely 

(b) For points (i+1) to (z): Averaging rows (i) and 
(iii); Dist=2*dely 

For these grids, though the adjacent pixel rows, both on top and bottom to the 
pixel row of interest are available, difficulty arises due to non equal number of pixels 
over these individual rows. This difference in the number of pixels is due to the 
variation in the size of the disc as traversal is done from top to bottom. Referring to Fig. 
3.11(b), it can observed that for a particular section of a grid, if a complete set of three 
pixels one over the other is not available, the averaging of associated parameters is to 
be performed over only two rows, which for the disc happen to be the central and the 
bottom rows. 

BOTTOM-MOST GRID : 

Grid structure: Row (iii) missing due to absence of Physical Geometry. 

Row (i) and (ii) available. 

Row (ii) 


Row of interest: 
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Averaging row/Dist.: Averaging rows (i) and (ii); Dist=dely 
The bottom-most grid is again constituted of only two rows, as is shown in Fig. 
3.1 1(c). Utilizing the symmetry of the model as only a quarter section is being analyzed, 
the physical geometry and hence the photoelastic data is not available below the line of 
symmetry. Due to this feature, the averaging of the associated parameters is done over 
the line of symmetry and a row of pixels just above it. 

3.8 Digital implementation of stress separation technique. 

As has also been mentioned earlier that the present analysis adopts theoretical data, 
partly realistic data and completely realistic data, sequentially to predict the nature of 
variation of stresses along x and y directions, for the whole field. Theoretical data is one 
which is obtained by using the standard formulations based on "theory of elasticity" 
(refer appendix I). The partly realistic data includes ( N) obtained from the experimental 
technique called phase shifting while completely realistic data refers to that combination 
where ( N) is obtained from phase shifting and (0) from either of the circular polariscope 
or the plane polariscope arrangement. 

3.8.1 Verification of algorithm, using theoretical data 

Consider any general pixel (P) in the region of interest as shown in Fig. 3.12. 

From the boundary information, it’s location in terms of pixels is known, say (px.py). 
The centre being taken as the origin, it’s spatial co-ordinates (x,y) can be generated as. 

x = px * delx and y — py * dely 



61 



Fig. 3.12 Generation of spatial coordinates for a pixel in the region of interest. 


With the spatial co-ordinates known for each pixel in the domain, corresponding 
fringe order ( N) and isoclinic parameter ( d ) can be obtained by using the standard 
formulations for disc under diametral compression (refer appendix I). On this basis, data 
for the entire field can be obtained and hence stress separation be performed. With the 
data available for the entire field, the next concern lies in determining the starting value 
for the integration procedure as proposed by shear difference method and Trebuna’s 
method. For this purpose the concept of stress tensor at the free boundary is applied. 

3.8. 1.1 TRANSFORMATION LAW 

For 2-D case, the stress tensor along x > -y > can be obtained by transforming the stress 
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tensor along x-y axis (Fig. 3.13). 




X 

y 

x’ 

cosd 

sin0 

y 

-sin0 

cosd 


Fig. 3.13 Transformation of stress along a new set of axes. 
Hence the transformation matrix [T] is 


m 


cos 6 sind 
-sin 6 cos 6 


(3.9) 


According to stress transformation law, the transformation matrix can be pre-multiplied 
while its transpose post-multiplied to the stress tensor along x-y direction to obtain the 
stress tensor along x'-y' direction. 


3.8.1. 2 THE FREE BOUNDARY 


When a boundary of the model is not loaded directly, it is called a free boundary. The 
normal and shear stresses on a plane tangential to a free boundary are therefore zero. 
The principal stress axes are normal and tangential to the boundary (Fig. 3.14). Hence 
one of the principal stresses, say a 2 , is zero while the non-vanishing principal stress a, 
is to be found out. 

As the direction of rotation here is reverse of the (+) convention, the sign of (6), 
needs to be reversed. Also a 2 and o n equals to zero. Hence the transformation law 


shapes up as: 
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cosd -sin 6 


°i 
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COS0 

sin0 

> °y. 


sin0 cos 8 


0 

0 


-sin d 

cosd 


yielding 

a x = o, cos 2 (6) 
a y = a, sin 2 (0) 

where: 

a, = N fjt ; a 2 being zero. 

N : fringe order on the boundary. 
f 0 : Material fringe value, 
t : thickness of the model 


(3.10) 



Fig. 3.14 The free boundary. 


It is realized by now that the stress tensor at the free boundary can be directly obtained 
from photoelastic data, if the isochromatic parameter (N) is known at the boundary. In 
this light, the shear difference method and Trebuna’s method are employed for analysis 
and the results verifying the algorithm, have been presented in Fig. 3.15 - Fig. 3.17. It 
can be seen that the results obtained by shear difference scheme and Trebuna s scheme 
are in close conformance with the theoretical results. 
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Fig. 3.15 


Use of theoretical data. Variation of stresses over horizontal scan line 
at y = 0 (mm), (a) variation of o x , a y (b) variation of t xy 






Fig 3 16 Use of theoretical data. Variation of stresses over horizontal scan une 
at y = R/2 (mm), (a) variation of ct x , a y . (b) variation oft,. 











4.6 


67 


Fig. 3. 



(b) 


(c) 




8 Variation of Isochromatic parameter (N), obtained f ror “ 
technique, for horizontal scan lines at (a) y - 0 (mm) (b) y ( ) 

(c) y = 3R/4 (mm) 
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3.8.2 Use of N from phase shifting; 6 from theory. 

Fringe order ( N) at every pixel can be obtained by digital processing of phase shifting 
images, as previously discussed. For an estimation of the accuracy of ( N ), so obtained, 
comparison with the theoretically generated (AO has been shown in Fig. 3.18, for three 
different horizontal scan lines. For the disc under diametral compression, the boundary 
with a known fringe order of zero can also be used as a reference for determining the 
accuracy of (AO from phase shifting technique. Figure 3.19 presents the value of (AO, 
obtained from phase shifting against theoretically generated (AO, along the boundary. 



Fig. 3.19 Variation of Isochromatic parameter (N), along the boundary of the 
quarter disc. 


Isoclinic parameter (0) can still be generated based on theoretical formulation for each 
pixel. This set of data can be employed to perform, a partially realistic analysis of stress 
field. The results presented in Fig. 3.20 - Fig. 3.22, for three different horizontal scan 
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Fig. 3.20 Use of N: phase shifting; 0: theoretical. Variation of stresses over horizontal 
scan line at y = 0 (mm), (a) variation of o x , cr y . (b) variation of 








Fig. 3.21 Use of N: phase shifting; 0: theoretical. Variation of stresses over horizontal 
scan line at y = R/2 (mm), (a) variation of a x , cr y . (b) variation of 
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Fig. 3.22 Use of N: phase shifting; 0: theoretical. Variation of stresses over horizontal 
scan line at y = 3R/4 (mm), (a) variation of a x , a y . (b) variation of 
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lines, reveal that though the level of accuracy in the stresses, separated by shear 
difference method and Trebuna’s method has fallen with the usage of phase shifting (TV), 
yet, the results are reasonably accurate. This comparison is of significant importance as 
it helps through, a step ahead in establishment of phase shifting technique as a reliable 
experimental technique, for determination of isochromatic parameter (jV). 

3.8.3 Use of N from phase shifting; 0 from plane polariscope. 

This approach focuses the use of completely realistic data for stress analysis. Here apart 
from obtaining the isochromatic parameter ( N) experimentally, the isoclinic parameter 
(0) is also obtained experimentally, using plane polariscope arrangement. The intensities 
of the recorded images are processed, as mentioned in previous sections and both the 
isochromatic parameter ( N) and the isoclinic parameter (0), obtained. As (0) is not 
defined for fringe orders of 0, 1,2 etc., excessive noise is present in the corresponding 
zones. A comparative analysis of (0) obtained from circular polariscope arrangement and 
the plane polariscope arrangement, has been presented for three different horizontal scan 
lines in Fig. 3.23. It can be observed that (0) obtained from the later is more accurate. 
It is for this reason that (6) from plane polariscope is further suggested for analysis. 

3.8.4 Use of N from phase shifting; 0 from plane polariscope; use of 
multiple loads. 

It has been observed that the accuracy of isoclinic determination is very poor near full 
fringe orders. In general, discontinuous isoclinic fringes are obtained. To bridge the gap. 
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Fig. 3.23 Comparative study of Isoclinic parameter 9: circular polariscopejs. plane 
polariscope arrangement, for various horizontal scan lines, (a) y — (mm) 
(b) y = R/2 (mm) (c) y = 3R/4 (mm). 
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Fig. 3.24 Comparative study of Isoclinic parameter 8, obtained from plane polariscope 
arrangement: Single load vs. Multiple load, for various horizontal scan lines, 
(a) y = 0 (mm) (b) y = R/2 (mm) (c) y = 3R/4 (mm) 
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use of interpolation approach 17 or multiple wavelengths 18 is reported. Extending the idea 
of multiple wavelengths, Petrucci 19 used white light to illuminate the model and 
extracted isoclinics using the red, green and blue intensity values of the colour image 
recorded by a colour image processing system. Then Ramesh and Mangal 15 proposed a 
new scheme, using two different loads to get continuous isoclinic fringes. Using this 
approach, first, the intensity information from the images corresponding to the load of 
335.05 N is processed and isoclinic parameter obtained. Then, the procedure is repeated 
for the raised load of 502.5 N. Finally, isoclinic for the field is obtained by perfor ming 
a logical ’OR’ operation between the two sets of isoclinic parameter. Figure 3.24 
presents a comparative study of isoclinic obtained using single load and that obtained 
using two loads, for various horizontal scan lines. Isoclinic obtained using two loads, 
is found to be largely improved over the single load isoclinic. Now with the fact 
established that the value of (0) can be best obtained from usage of multiple loads, with 
plane polariscope arrangement, it’s usage for further analysis is suggested, employing 
the same technique. 

Figure 3.25 shows the variation of (0), obtained from use of multiple loads, along 
the boundary of the quarter of the disc. It can be observed that the value of 0, along the 
boundary is a constant whereas the actual variation of 0, along the boundary, should be 
of the form, as represented by the theoretically generated 0. This difference can be 
attributed to the fact that 0, given by the eq. (3.3) or eq. (3.5), may be indeterminate, 
at some points. While processing the intensity information of the grabbed images and 
obtaining the values of isoclinic parameter 0, for the whole field, the value at such 
indeterminate points, was arbitarily assigned. The constant value of 0, seen along the 
boundary is nothing but this arbitarily assigned value itself. 
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Fig. 3.26 Use of N: phase shifting; 0: plane polariscope / multiple loads. Variation 
of stresses over horizontal scan line at y = 0 (mm), (a) variation of a x , 
cr y . (b) variation of x xy 
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(*) 



(b) 


Fig. 3.27 Use of N: phase shifting; 0: plane polariscope / multiple loads. Variation 

of stresses over horizontal scan line at y = R/2 (mm), (a) variation of a x , 
cr r (b) variation of T xy 
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Fig. 3.28 


Use of N: phase shifting; 0: plane polariscope / multiple loads. Variation 
of stresses over horizontal scan line at y = 3R/4 (mm), (a) variation of 
a x , a y . (b)variation of x xy 
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(b) 




1 

Fig. 3.30 Variation of photoelastic data for horizontal scan line at y = R/4 (mm). 

(a)variation of N: phase shifting, (b) variation of 0: plane polanscope 

arrangement using multiple loads. 
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The results for various horizontal scan lines have been presented in Fig. 3.26 to Fig 
3.29. It is seen that the stresses evaluated using completely realistic data are not in 
agreement with the theoretically evaluated stresses. However, the results for the 
horizontal scan line at y = R/4, are close to the theoretical results. This can be accounted 
to the better input data, for this particular scan line. The photoelastic data in use for this 
line is presented in Fig. 3.30. 

It is important to realize that the indeterminacy of 0 is not restricted to the 
boundary alone. It could be anywhere across the field also but the indeterminacy of 0, 
at the boundary has relatively more effect in terms of reducing the accuracy of final 
solution. If, the value of 0 at the pixel on the boundary which also happens to be the 
starting point for the integration procedure as proposed in shear difference scheme or 
in Trebuna’s scheme, is improved then the corresponding solution can be expected to 
be better. One possible improvement could be in the form of assignment of a well 
defined value of 0, closest to the boundary, as the value on the boundary itself. 
Referring to Fig. 3,1 1(a), assume a general case where 0 at the pixel boundary, say P, 
is indeterminate and the closest pixel having a well defined value of 0 is Pj. If the value 
of 0 at P; is assigned to all the pixels towards the boundary, the corresponding solution 

i 

should be better. Figure 3.31(b), shows the modified 0, along the boundary. The 
corresponding solutions for arbitrary horizontal scan lines are presented in Fig. 3.32 
It can be observed that these results are relatively better than those corresponding to the 
un-modified 0. Also, the results corresponding to Trebuna’s method seem to have 
improved more than those corresponding to shear difference method. 



Distance along y axis (mm) 



Fig. 3.31 


Modification of Isoclinic parameter ( 6 ), along the boundary (a) the 
proposed scheme (b) improved variation along tire boundary. 
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3.8.5 Whole field representation of the results. 


All the results presented thus far have been for a particular horizontal or a vertical scan 
line. Now, the aim is to represent the variation of a x , a y and x xy , for the whole field. The 
eq. 2.1 relates the difference of principal stresses at a point to the fringe order at that 
point. Analogous to this, if (a, - a 2 ) is replaced by cr x , then (N) obtained from the 
equation would correspond to <j x . If ( N) corresponding to a x is represented as N a , then 
it can be given by the relation: 



(3.11) 


If N 0 is plotted, it would be in one sense the field representation of ct x itself, the 
constant (t/F a ) being the only discriminatory factor between the two. Similarly the plots 
for N a and N x can also be obtained. Fig. 3.33(a), Fig. 3.33(b) and Fig. 3.33(c) represent 
the variation of a x , a y and r xy respectively. This solution set is obtained by utilizing 
theoretically generated photoelastic data, with the standard formulations for the disc 
under diametral compression. The main emphasis of the thesis apart from the digital 
aspects of stress separation schemes has been to investigate the trend with which the 
accuracy of the results vary, as the input data is changed from theoretical to partly 
experimental and then to completely experimental. Figure 3.34 shows this trend for ct x . 
Figure 3.34(a) and Fig. 3.34(b) represent the results obtained from shear difference 
method and Trebuna’s method respectively, using the theoretically generated data. 
Figure 3.34(c) and Fig. 3.34(d) show the results corresponding to the combination of 
phase shifting ( N) and theoretical ( 0 ). The results can be observed to have detiorated. 
Finally, the resort being made to completely realistic data, wherein (0) is also 



cr x : Theory 
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CT y : Theory 



(b) 


Txy : Theory 



Fig. 3.33 Pseudo fringe patterns showing variation of stresses, determined using 
theoretical data and standard solution for disc under diametral 
compression. 
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FIg - 3 ' 34 00 ^ftecry; « - Trebuna 

e N p-sftr’« theo^- shear diff. (d) N:p-sftg; (htheory - Trebuna 
t) N'p-sftg; S:mult.loads - shear diff.(f) N:p-sftg; (hmulUoads-Trebuna 



Shear difference 


Trebuna 




Shear difference 


Trebuna 




m 


Shear difference 






Trebuna 


m2 


Pseudo fringe patterns for cr y for various sets of input data. 

(a) N:theory; 0: theory - shear diff. (b) N:theory; 0:theory - Trebuna 
(c) N:p-sftg; fhtheory - shear diff. (d) N:p-sftg; 0:theory - Trebuna 
(f) N:p-sftg; 0:mult.loads - shear diff.(f) N:p-sftg; 0 :mult. loads-Trebuna 








Shear difference 


Trebuna 
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(a) 


Shear difference 



(c) 


Shear difference 



(e) 



(b) 

Trebuna 



(d) 

Trebuna 


(f) 



Fig. 3.36 


Pseudo fringe patterns for t xy a JJ° u ^ S J tS ^ f 0 P theory - Trebuna 
(a) N '.theory; theory - she*- *ff. Trebuna 

(c) N:p-sftg; 0 '.theory - shear diff. () -P | 0 - mu lt.loads-Trebuna 
(f) N:p-sftg; (hmult.loads - shear diff.(f) N.p sttg, b.n 
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experimentally obtained , using multiple loads with plane polariscope arrangement, the 
corresponding results have been presented in Fig. 3.34(e) and Fig. 3.34(f). These results 
reflect the inaccurate and noisy (6). Figure 3.35 and Fig. 3.36 present the variation of 
CT y and x x y, respectively. It is worth mentioning here that this mode of fringe plotting is 
not a true representation of the results, due to the limitation of the plotting code. The 
stresses being plotted as fringes are only those which lie within a user defined tolerance 
zone. Any stress value even slightly off this zone, is not plotted, giving a haphazard and 
a noisy view. Despite of this fact, the fringe representation of the stresses is employed 
to give a feel of stress variation over the entire field. 

3.9 Plate with a hole 

The previous sections have discussed the establishment of shear difference scheme and 
Trebuna’s scheme as reliable schemes for stress separation. The trend with which the 
accuracy of the results vary, as the input data is sequentially changed, first from 
completely theoretical to partly experimental and finally to completely experimental, has 
also been discussed. This section discusses the application of above mentioned stress 
separation schemes over a plate with a hole, subjected to uniaxial tension. The data 
employed for analysis is completely experimental where (N) is obtained from phase 
shifting and (0) from plane polariscope arrangement. 

3.9.1 Boundary Information 


It has earlier been mentioned that in order to detect the boundary of the region of 



91 


interest, the maximum grey level over the entire image is detected. The variation of 
these grey levels can be from 0 to 255. Owing to the fact that the maximum grey level 
over the image would be less than or equal to 255, any pixel with grey level value of 
255 would represent the boundary pixel. Adopting this principle, to mark the boundary, 
Sreedhar[J had developed some primitives. These are line, circle, ellipse and 3 -point arc. 
Utilising the symmetry of the plate with the a hole, only a quarter section needs to be 
evaluated. It is to be realised that the complete boundary of the quarter section of the 
plate, under consideration, cannot be directly determined using any of the primitives 
mentioned above. For determination of the complete boundary, first the hole in the plate 
is considered. The boundary of this hole can be detected by simulating it as an ellipse 
and using the primitive developed by Sreedhar. The spatial location of the three 
remaining comers of the plate can be known from the grabbed image, in terms of pixels. 
This information in conjunction with the known hole boundary can be used to determine 
the complete boundary of the plate with a hole. 

As it has been proved in earlier sections that the photoelastic data obtained from plane 
polariscope arrangement is superior to that obtained from circular polariscope 
arrangement, images of the plate under consideration (refer appendix III) are grabbed 
using the plane polariscope arrangement. The sequence of six images grabbed are shown 
in Fig 3.37. The details regarding the determination of isochromatic parameter (N) and 
isoclinic parameter (6) have bben discussed in section 3.3.2. Now, with the boundary 
detected and the experimental data available for the entire field, the concern lies in 
determining the starting point for initiation of integration procedure as proposed in shear 
difference and Trebuna’s scheme. For this purpose, again the free boundary is made use 
of. A boundary which is not directly loaded, is referred as a free boundary. The normal 
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Fig. 3.37 The sequence of six images of a plate with a hole under uni-axial tension 
recorded for phase shifting technique, by plane polariscope arrangement. 
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and shear stresses on a plane tangential to free boundary are therefore zero. Also, as the 
principal stress axes are normal and tangential to the boundary, one of the principal 
stress is zero on the free boundary. Finally, (refer section 3. 8. 1.2) ct x and cr y on the free 
boundary are given by : 

a x = cr, cos 2 (0) 

Oy = ct, sin 2 (0) 

CTj = N fji ; ct 2 being zero. 

N : fringe order on the boundary. 
f c : Material fringe value, 
t : thickness of the model 

Now, the stress separation schemes are applied and results obtained. Ramesh, Mangal 
and Pathak 20 , had solved the problem of the plate with a hole, using finite element 
analysis and presented the variation of principal stresses a, and ct 2 in the form of fringes. 
So, as to present a comparative analysis, the stresses separated by shear difference 
scheme and Trebuna’s scheme have been converted to the principal stresses and plotted 
in the form of the fringes (refer section 3.8.4) and presented in Fig. 3.38. Figure 3.38(a) 
and 3.38(b) show ct, and ct 2 contours, obtained from finite element analysis. The fringe 
contours for the principal stresses obtained by shear difference scheme are shown in Fig. 
3.38(c) and 3.38(d) while those obtained by Trebuna’s scheme are presented in Fig. 
3.38(e) and 3.38(f). With, erroneous (6), Trebuna’s method is expected to offer better 
results, due to the averaging effect of (0), across the grid lines. It can be clearly seen 
from the graphs that shear difference’s solution has certainly been improved by 
incorporating Trebuna’s modification to it. 
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form of pseudo fringe patterns. 

From the results, for the disc under diametral compression, it is revealed that the 
results obtained by Shear difference method and Trebuna’s method by employing 
theoretical data are in complete agreement with those obtained by employing the 
standard formulations, based on theory of elasticity. The results obtained corresponding 
to usage of phase-shifting (N) and theoretical 0, are reasonably accurate. This establishes 
phase-shifting as a reliable technique for determining the isochromatic parameter (N). 
The indeterminacy of (0), at full fringe orders, poses problems and these points are 
reflected as noise points. The usage of erroneous (0) along with the phase-shifting (N) 
leads to extremely poor results which are in no way reliable. It is shown that the 
averaging effect of (0) in Trebuna’s method causes improvement in results over shear 
differences’s results. This effect is also reflected in the stress contours for the plate with 
a hole, subjected to uniaxial tension. 

4.2 Suggestions For Future Work 

1. Since the role of 6 plays an important role on the overall accuracy , 
newer technique for 6 evaluation by phase shifting has to be looked into. 

2. Recent studies have indicated that the evaluation of N by phase shifting 
tec hni que can be improved further. This methodology needs to be 
incorporated for stress separation. 

3. The method should be applied for a variety of problems to improve the 
ruggedness of the solution. 
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APPENDIX-I 

THEORETICAL STRESS FIELD FOR DISC 
UNDER DIAMETRAL COMPRESSION 


The stress field in a circular disk under diametral compression can be obtained 
using theory of elasticity from Bousinesque’s solution for a compressive load on an half 
plane and then applying the principle of superposition appropriately. The stress field is 
obtained as 


where. 




(R-y)x 2 , (R+y)x 2 1 1 

f 


V 

* — 

°x 

- JE, 

T it 

(R-y) 3 . (R + y) 3 1 

°y 

r 4 r 4 D 1 

r l r 2 




i *y j 


(R+y) 2 x (R-yfx 



4 4 

H r i 


r 

\=X 2 +{R-j) 2 


rl=x 2 +(R+yf 


P: is the applied load (N). 

R: is the disc radius (mm). 

D: is the disc diameter (mm), 
t: is the thickness of the disc (mm). 
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APPENDIX-U 

MODEL : DISC UNDER DIAMETRAL 
COMPRESSION 


PHYSICAL DIMENSIONS: 


(i) Diameter : 

(ii) Thickness : 

60.06 mm 

5.04 mm 

MATERIAL PROPERTY : 


Material fringe value F CT : 

1 1.64 N/mm/fringe 

IMAGE FRAME SPECIFICATIONS : 


Image used 

: dskpltl2.p01 

: dskpltl2.p02 

: dskpltl2.p03 

Centre 

: (127,127) 

Increment along X-axis, delx 
Increment along Y-axis, dely 

: 0.3532941 mm 

: 0.2461475 mm 

Aspect ratio 

: 0.6967213 
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appendix-iii 

MODEL : PLATE WITH A HOLE UNDER UNI- 
AXIAL TENSION 


PHYSICAL DIMENSIONS: 


(i) Width 

(ii) Length 

(iii) Thickness 

(iv) Hole Diameter: 


37 mm 
174 mm 
5.26 mm 
60.06 mm 


MATERIAL PROPERTY : 

Material fringe value F 0 


11.5 N/mm/fringe 


IMAGE FRAME SPECIFICATIONS . 

Image used 


platewl5.p01 

platewl5.p02 

platewl5.p03 


Centre 


(154,225) 


Increment along X-axis, delx 
Increment along Y-axis, dely 


0.362745 mm 
0.257549 mm 


Aspect ratio 


0.710 


